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Abstract. The now extinct ^°^Pb is a pure s- process ra- 
dionucHde (^1/2 = 1-5 10'' y) of possible substantial cos- 
mochemical interest. As a necessary complement to the 
detailed theoretical study of the nuclear physics and as- 
trophysics aspects of the ^"^Pb - ^os^p]^ pg^jj. QQ^j-ried out 
by Yokoi et al. (1985), and to the recent calculation of 
the 205p]3 production in Wolf-Rayet stars by Arnould et 
al. (1997), this paper addresses for the first time in some 
detail the question of the survival of this radionuclide in 
thermally pulsing AGB stars. This problem is made diffi- 
cult by the high sensitivity to temperature and density of 
the rates of the weak interaction processes that are able 
to produce or destroy ^°^Pb. In view of this sensitivity, a 
recourse to detailed stellar models is mandatory. With the 
help of some simplifying assumptions concerning in par- 
ticular the third dredge-up characteristics, some of which 
(like its depth) being considered as free parameters, pre- 
dictions are made for the ^°^Pb contamination of the stel- 
lar surface at the end of a pulse-interpulse cycle following 
a series of a dozen of pulses in three different intermediate- 
mass stars (M = 3M0,Z = 0.02; M = 6Mq,Z = 0.02; 
M = 3Mq,Z= 0.001). It is concluded that the chances 
for a significant ^"^Pb surface enrichment are likely to in- 
crease with M for a given Z, or to increase with decreasing 
Z for a given M. More specifically, following the consid- 
ered pulses at least, the enrichment appears to be rather 
unlikely in the 3Mq star with Z = 0.02, while it seems to 
be much more probable in the other two considered stars. 
It is also speculated that the {3Mq,Z = 0.02) star could 
possibly experience some ^°^Pb enrichment following later 
pulses than the ones considered in this paper. 
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1. Introduction 

There is at present ample observational evidence that 
some meteoritic material of solar system origin carries the 
signature of the in-situ decay of by now extinct radionu- 
clides with half-lives from about 10^ to approximately 10^ 
y. These observations have far-reaching implications for 
our knowledge of the history of the forming solar sys- 
tem, and provide in particular a key information about 
the time span between the last nucleosynthetic events that 
have modified the composition of the solar nebula and the 
formation of the solar system solid bodies. For a recent 
review, the reader is referred to e.g. Swindle (1993) and 
references therein. 

Among the extinct radionuclides of potential cosmo- 
chemical interest, ^°^Pb {ti/2 = 1.5 lO'' y) exhibits the 
distinctive feature of being of pure s-process nature. If its 
signature in meteorites could be identified, it would thus 
very usefully complement the information provided by the 
other extinct radionuclides. Unfortunately, its abundance 
at the start of the solidification process in the early solar 
system is very poorly known. In fact, just an upper limit 
on this quantity is available (Huey & Kohman 1972). 

A detailed theoretical study of the nuclear physics and 
astrophysics aspects of the ^°^Pb - ^o^Tl pair has been car- 
ried out by Yokoi et al. (1985, hereafter YTA85). These 
authors have shown that, in contradiction with an earlier 
estimate (Blake & Schramm 1979), enough 205pj-, could 
be produced in certain s-process conditions to allow the 
development of a s-process chronometry. This of course re- 
quires a renewed search for the signature of this radionu- 
clide in meteorites (a step in this direction has been taken 
by Chen & Wasserburg 1994), as well as a reliable enough 
prediction of the ^"^Pb yields from all possible s-process 
sites. 

In order to explore the latter question, YTA85 have 
considered simple parametric models that can at best 
mimic the complex conditions prevailing in the stellar in- 
teriors. The first quantitative attempt to evaluate ^"^Pb 
stellar yields has been conducted by Arnould & Prantzos 
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(1986) on grounds of detailed models of massive mass- 
losing stars of the Wolf-Rayet type. Their early predic- 
tions, confirmed recently by Arnould et al. (1997), re- 
inforce the view already expressed by YTA85 that in- 
deed ^"^Pb could be a valuable s-process chronometer. 
Asymptotic Giant Branch (AGB) stars, on the other hand, 
have also been considered by YTA85 as potential sites for 
the ^"^^Pb enrichment of the interstellar medium (see also 
Wasserburg et al. 1994). Their estimate, however, is based 
on a very schematic representation of the recurrent ther- 
mal pulse - interpulse cycles characterising these objects. 
In addition, the very important question of the survival 
of 205p|-, during the interpulse phases has been left unan- 
swered, nor has it been tackled by Wasserburg et al. (1994) 
either. 

The aim of this paper is to fill this gap by following 
quantitatively the decay of ^°^Pb in essentially neutron- 
free locations of AGB stars during their interpulse phases. 
These computations are based on the /3-transition rates 
derived by YTA85 and on detailed AGB models for three 
different stars (M = SM©, Z = 0.02; M = GM©, Z = 
0.02; and M —3 M©, Z = 0.001) whose layers experiencing 
the pulse - interpulse episodes are supposed to be enriched 
with s-process products, and in particular with ^"^^Pb. It is 
then assumed that the third dredge-up (hereafter 3DUP) 
is responsible for the transport of the surviving ^"^Pb to 
the surface of the AGB stars before being ejected into the 
interstellar medium by their winds. 

Some generalities about the AGB features of interest 
here (pulses and s-process nucleosynthesis) are briefly re- 
viewed in Sect. ||. The specific aspects of the /3-decay of 
205p^ in AGB conditions are described in Sect. |[ and 
the prescriptions adopted for the calculation of its abun- 
dance evolution in neutron-free AGB layers are presented 
in Sect. ^ The predictions of the survival of ^ospj^ 
the three considered model stars are presented in Sect, [s] 
and discussed in Sect. ^. Some conclusions are drawn in 
Sect. 0. 

2. The AGB stars 
2.1. Thermal instabilities 

The structural evolution of the H- and He-burning shells 
(denoted hereafter HBS and HeBS, respectively) during a 
thermal instability (pulse) characteristic of the AGB stars 
is displayed in Fig. |l|. For the sake of further discussion, 
four distinct cycle phases are identified, the main charac- 
teristics of which are summarized below (more details can 
be found in Mowlavi 1995): 

Phase 1: The interpulse phase. Hydrogen burns in a thin 
layer of a few 10"'' M© located at mass fractions in- 
creasing with time. The H-burning products, including ^^C 
and ^^N, are left behind in the intershell region (layers be- 
tween the HBS and the HeBS). The HeBS, whose mass is 
about 10^^ M©, is almost extinct, and the luminosity of 



the star is essentially provided by the HBS. As the mass of 
the He-rich layers increases, a thermal instability is trig- 
gered in the HeBS (Schwarzschild & Harm 1958), leading 
to the development of a convective pulse; 

Phase 2: The growth of the convective pulse. A convective 
zone develops just above the location of the maximum en- 
ergy production of the HeBS, its upper boundary coming 
close to the H-rich layers. The ashes left behind by the 
HBS during the interpulse period are engulfed into the 
pulse. The pulse duration is about 500 times shorter than 
the interpulse one; 

Phase 3: The decay of the pulse. The energy accumulated 
in the HeBS escapes through the now extinct HBS, and 
the convective pulse recedes; 

Phase 4: The afterpulse phase. Following the pulse, the 
structural and energetic evolution is complex and domi- 
nated by thermal relaxations. The temperature and den- 
sity in the He-rich layers first drop as a result of the expan- 
sion of the intershell region, and then increase due to the 
structural readjustment. At the same time, the convective 
envelope reacts to the thermal relaxations by deepening 
into the H-depleted layers, possibly reaching zones con- 
taining the ashes of the convective pulse. The mixing of 
these ashes into the envelope and their transport to the 
surface constitute the 3DUP. This dredge-up is required 
by the observation of the surface composition peculiari- 
ties of AGB stars, and is predicted by some models after 
a sufficient number of pulse/interpulse cycles (Mowlavi 
1997). The model pulses dealt with in this paper (Sect. 
5) are not accompanied with the 3DUP, which appears to 
develop only at later pulses in the considered stars. How- 
ever, some observations suggest that the 3DUP might oc- 
cur much earlier than predicted (Mowlavi et al. 1996). In 
view of the uncertainties still affecting the precise time 
of occurrence and extent of the 3DUP, we simply assume 
in the following that it indeed takes place at the time of 
deepening of the convective envelope in all the pulses to 
be discussed later. This time also defines the end of the 
afterpulse phase, and correspondingly the start of the fol- 
lowing interpulse period. 

2.2. The s-process 

While observation demonstrates that the s-process can 
develop in AGB stars, its modelling remains a warmly 
debated question. This concerns more specifically the 
precise mechanism of production and the amount of 
neutrons made available for the process. The two clas- 
sically envisioned neutron sources, ^'^C {a , n) ^^O and 
^^Ne (a , n) ^^Mg , have an efficiency that is expected to 
depend very much on the considered stars, as well as on 
the precise evolutionary phase of a given star. The relia- 
bility of the predictions in this field is drastically limited 
by astrophysics and nuclear physics uncertainties (Drotleff 
et al. 1993, Frost & Lattanzio 1995). 
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Fig. 1. Structural evolution of the intershell regions of the 3Mq star with Z — 0.02 around the 12th pulse. The origin of the 
abscissa is arbitrarily set to time to, the value of which is indicated in the figure. Filled regions correspond to convective zones. 
The dotted lines identify the location of maximum energy production in the H-burning (top) and He-burning (bottom) layers. 
The different phases of an instability cycle described in the main text are also indicated 



In the models considered in the present paper, the tem- 
peratures limit the neutron production to ^'^C (a , n) ^^O . 
The available ^^C is found to come from the CN cycle ex- 
clusively. This amount is well known to be insufficient for 
producing the neutrons required for the development of a 
full s-process. This problem classically encountered in the 
modelling of the s-process in AGB stars is circumvented 
by assuming that some protons from the envelope are 
brought "semi-convectively" into the underlying ^^C-rich 
layers during the afterpulse phase (Iben & Renzini 1982). 
Recent evolutionary calculations (Herwig et al. 1997) pre- 
dict similar mixing to occur by convective overshoot at 
the base of the envelope. In both cases, after reignition 
of the HBS in the early interpulse phase, some ^^C could 
be produced by (p , 7) ^^C . Part at least of this 
of "primary" nature could be subjected to ^^C (a , n) ^^O 
either during the interpulse phase or later after ingestion 
in the subsequent pulse. These two burning modes will 
be referred to in the following as the "Interpulse Scenario 
(IS)" and the "C Pulse Scenario (CPS)" , respectively. Of 
course, the s-process abundances emerging from a pulse 
may result from the complementary and sequential action 
of the IS and of the CPS. We define this scenario as the 
IS-hCPS. 

Even if the ^^C (a , n) ^^O neutron source is currently 
believed to be responsible for the s-process in AGB stars, it 
remains that the amount of neutrons that can be produced 
in such a way is still highly uncertain, and its dependence 
on stellar mass and metallicity is unknown. 



3. The production and destruction of ^osp]^ 

As emphasized in Sect. 1, ^''^Pb has the distinctive fea- 
ture of being of pure s-process origin, and could thus be 
synthesized in the IS and/or CPS. The produced quanti- 
ties depend of course on the amount of available neutrons, 
which cannot be estimated self-consistently and reliably at 
the present stage of development of the models. 

A complementary question, which has not been ad- 
dressed properly yet, concerns the survival of ^"^^Pb in the 
essentially neutron-free layers in which the radionuclide 
may reside after its production by the s-process. Obvi- 
ously, the level of ^°^Pb destruction in such conditions 
depends on the relative values of the relevant stellar evo- 
lutionary lifetimes and of the timescales of production or 
destruction of this radionuclide through weak interaction 
processestl. As demonstrated by YTA85, the latter hfe- 
times may depend drastically on the stellar conditions. 
This makes the question of the survival of ^osp]^ non- 
trivial problem. More specifically, in a stellar plasma the 
thermal population of low-lying nuclear excited states, 
ionization, and the possibility of capture of free electrons 
can modify substantially the experimentally known ^"^Pb 
half-life, and even open the possibility for the laboratory 
stable ^os^pj become /3-unstable and to transform into 
205pj^^ Figure || illustrates the impact of the stellar con- 



^ In the considered neutron-free stellar environments, 
charged particle reactions or photodisintegrations cannot af- 
fect the ^°^Pb abundances. 
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Fig. 2. YTA85 rates Ae of ^"'^Pb electron capture and A_ of 
20Brpj ^.jjggg^y versus temperature for the electron number den- 
sities iVe (in cm"^) of lO^'^ (full Hnc), lO^*^ (dot-dashed line) and 
10^^ (dashed line). The uncertainties of nuclear origin affect- 
ing these predictions are discussed in detail by YTA85. They 
concern mainly some experimentally unknown /3-decay matrix 
elements, as well as the ^"^Pb - ^°^T1 mass difference 



ditions on the ^°^Pb and ^os-pj decay rates. Let us just 
emphasize that 

(1) the ^°^Pb e~-capture half-hfe gets much shorter than 
its terrestrial value of 1.5 10'' y at temperatures that are 
high enough (typically in excess of a few 10^ K) for the 
first excited state of ^°^Pb (excitation energy of only 2.3 

keV) to be significantly populated; 

(2) when its level of ionization is high enough, ^os^pj ^g^^^ 
transform into ^°^Pb via bound-state /3-decay, its corre- 
sponding half-life depending very sensitively on both tem- 
perature and density. Its lifetime is shorter than the ^°^Pb 
one for temperatures T ^ 1.3 10* K at a density p « 300 
gcm~^, this critical temperature increasing with increas- 
ing densities. 

On such grounds, YTA85 predict that ^°^Pb is in dan- 
ger of being transformed into ^os-pj electron captures 
in neutron-free regions at 10^ ^ T <, 10® K, the level 
of destruction being highly sensitive in this temperature 
range to the free electron number density and to the evolu- 
tionary timescales. In hotter locations, this destruction is 
prevented by the reverse transformation ^ostpj f^p- -j 205p|j _ 
Finally, at T ^ 10® K, the ^osp]-, e~-capturc lifetime gets 
longer than the typical intcrpulse period ( 10® y) ex- 
pected in AGB stars, so that 205p^ j^g^g time to decay 
in significant amounts before the occurrence of the next 
pulse. Of course, these temperature limits are only indica- 
tive since the rates are density dependent. 

In view of the above considerations, it appears hope- 
less to evaluate in a reliable way the level of destruction of 
205p|j in stellar neutron-free locations following the CPS 
or IS without recourse to a self- consistent abundance cal- 



culation conducted in the framework of AGB models. The 
first computation of this kind is described below. 

4. The stellar models and numerical assumptions 

The time variations of the ^°^Pb abundance is fol- 
lowed during some characteristic pulse - interpulse cy- 
cles through detailed model calculations of three AGB 
stars with different masses and metallicities referred to 
in the following as M3Z02 (M=3 M0,Z=O.O2), M6Z02 
(M=6M0,Z=O.O2), and M3Z001 (M=3 M0,Z=O.OO1), 
the evolution of these stars being computed all the way 
from the pre-main sequence phase. Details about these 
evolutionary models , as well as about the computer code 
being used can be found in Mowlavi (1995) and Mowlavi 
et al. (1996). 

These detailed AGB models fail to predict self- 
consistently the s-processing and the 3DUP, so that some 
assumptions have to be made in order to compute the 
205pj-, suj-face enrichment that may follow from a thermal 
pulse phase. More specifically, we suppose that 

(i) all the considered stars can produce ^°^Pb. This pos- 
sibility is substantiated by the fact that the considered 
model stars are found to produce neutrons through the 
CPS or IS-FCPS (see Sect. 5); 

(ii) in line with assumption (i), a ^'^"''Pb number density 
equal to unity is imposed in the whole convective pulse 
region at the moment of its maximum extension in mass, 
considered as the initial time tin in the abundance calcu- 
lations reported below. This choice of initial abundances 
is justified by the fact that we are not interested here in 
the absolute level of ^°'''Pb production, but rather in the 
extent of its survival after its synthesis; 

(iii) the 3DUP occurs at the time tdredge of deepest extent 
of the convective envelope (see Sect. 2.1). It is considered 
to mix homogeneously part of the remaining ^°^Pb into 
the whole convective envelope up to the surface. The frac- 
tion of ^"^^Pb that is effectively transported to the surface 
depends on the depth of the 3DUP, taken here as a free pa- 
rameter. In addition, the mixing is assumed to take place 
during a time interval that is short compared to the life- 
times of the ^°^Pb production or destruction mechanisms. 
These processes are thus frozen during the 3DUP; 

(iv) after the 3DUP, the ^°^Pb abundance is followed in the 
convective envelope over the time tinter of duration of the 
interpulse. In absence of a self-consistent 3DUP model, the 
physical conditions in the envelope are assumed to recover 
their pre-dredgc-up values; 

(v) in the tin < t < tdredge time interval, the layers 
loaded with ^°'''Pb following assumptions (i) and (ii) are 
neutron-free, so that the ^°^Pb abundance is governed 
by the simultaneous action of just ^°^Pb(e~z/)^°^Tl and 
2"5Tl(/3-)2n5pb. The phase of decay of the modelled 
pulses (see Sect. 5) is indeed found to be neutron free, 
the ^^C possibly engulfed in the pulse convective tongues 
being fully destroyed at f = tin, while ^^Ne cannot burn. 



N. Mowlavi et al: The survival of Pb in intermediate-mass AGB stars 



5 



On grounds of the above assumptions, the calculation 
of the evolution of the ^°^Pb abundance profile from tin 
up to the end of a given interpulse is then fully coupled 
to the detailed modelling of the pulse decay, afterpulse, 
and interpulse phases. For each stellar model, we limit 
ourselves to the consideration of a single pulse-afterpulse- 
interpulse cycle identified in Table 1 . This simplification is 
justified by the absence of self-consistent 3DUP phases in 
the computed models. It will be dropped in future models 
if they succeed to predict such a dredge-up. Even if the 
number of studied cycles is limited, they cover a quite large 
variety of situations ^°^Pb can face after its production. 

In radiative neutron-free layers, the ^''^Pb abundances 
follow from 

iV2a5pb(t + At) = (1) 
(^iV205pb(t) - 7V205^) exp'^^* +iV205^, 

with A = A_ + Ae, and 

^205 = A^205pb(t) + A^205Tl(t) 

= A^205pb(t + A<) + Ar205Ti(t At). (2) 

In these equations, Ni(t) is the abundance of nucleus i 
at time t, At the time step over which the abundance 
change is calculated. Following assumption (ii) of Sect. 4, 
N205piy{tin) = 1 and iV205xi(tin) = are adopted. On the 
other hand, A_ and Ag are the values of the ^os^j /3-decay 
and ^°^Pb e~'-capture rates appropriate for the temper- 
ature, density and composition (and thus electron con- 
centration) of the considered radiative layer. Equation (2) 
expresses the constancy of the total amount oi A = 205 
isobars in absence of any other transformation than the 
considered /?-transitions. 

In convective regions, instantaneous mixing is as- 
sumed. The evolution of their ^"^Pb content is still given 
by Eqs. (1) and (2) in which A_ and Ag are interpreted as 
mass averages of the local /3-decay and e "-capture rates 
over the whole convective region. 



Table 1. Some characteristics at the given pulse number of 
the three cases for which the ^"^Pb abundance calculations are 
performed. T{^^'^ is the maximum temperature reached at the 
base of the given pulse and Tbe the temperature at the base of 
the envelope during the following interpulse. The temperatures 
are expressed in units of 10^ K. The last column identifies the 
scenario assumed to be responsible for the ^°^Pb abundance 
production (see main text) 
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Fig. 3. (a) Structural evolution of the intershell regions dur- 
ing the decay of pulse 12 and during the afterpulse of the 
M3Z02 star. Filled areas have the same meaning as in Fig. ^ 
Long-dashed lines represent locations of constant temperature, 
labelled in units of 10^ K. Short-dashed lines indicate locations 
of constant density (in gcm~^). Values in parentheses indicate 
logarithms of the density. The origin of time is set at to, the 
value of which is displayed on the abscissa, (b) Same as (a) but 
for contours of equal ^"^Pb abundance. Numbers in parentheses 
refer to logarithms of the abundances 

5. Model predictions 

5.1. The M3Z02 case 

Up to at least pulse 12, CPS dominates the ^"'^Pb pro- 
duction. During the decay phase of this pulse, the ther- 
modynamic conditions in the convective tongue are such 
that ^"^Pb essentially survives. Though the typical (mass- 
averaged) lifetime of ^"^Pb is of the order of the « 30 y 



6 



N. Mowlavi et al: The survival of Pb in intermediate-mass AGB stars 



200 
150 
100 
50 

-8 
-9 
-10 
-11 



(a) 

3 Mq, Z=0.02 
after pulse 12 



log(p) 



I I I I I I I I I I I I I I I I 



■ I I T 



log(Xj 



^ X log(X_) 




\','}- I", I,', 'I,', V,'l', 'J',' 



>1 0.615 0.62 0.625 0.63 0.635 



Fig. 4. Mass profiles of some quantities in the M3ZQ2 star 
at the end of the 12th afterpulse, which occurs at t = 408 y 
in the time abscissa of Fig. |; (a) Temperatures (in units of 
10^ K) and densities (in gcm~'^), to be read on the left- and 
right-hand ordinates, respectively; (h) ^°^Pb e~-capture (A^) 
and ^"^Tl ^"-decay (A_) rates, expressed in s~^; (c) ^"^Pb and 
^"^Tl abundances. The helium mass fraction profile (dotted 
line) helps locating the envelope. The extent of pulse 12 is also 
shown 



duration of the decay phase, the lifetime of is k, 12 
y preventing great depletions of ■^°^Pb. 

In the afterpulse phase. Figs, ^('bj and ^(c) demon- 
strate that the destruction of ^°^Pb is severe in layers ex- 
periencing temperatures T <, 10^ K, as already predicted 
by YTA85 (see also Sect. This corresponds to about 
the upper half of the region formerly covered by pulse 
12. In this mass region, Ag is of the order of 10~^ s~^ [see 
Fig. §('bj], corresponding to a ^°^Pb lifetime of about 30 y. 
This is roughly ten times shorter than the afterpulse dura- 
tion, so that ^°^Pb decays largely into ^°^T1. On the other 
hand, the reverse transformation is too slow to impede 
this destruction. The situation in this respect is quite dif- 
ferent in deeper layers vacated by the pulse, where A_ can 
be comparable to Ag. In such conditions, the destruction 
of ^"^Pb is largely avoided, its abundance having time to 
reach its local equilibrium value iVa'o'iipj^ — 7V205A--/A. 

The results just described imply that a very deep 
3DUP of some 10~^Mq is required at the end of the af- 



0.9495 



0.949 



0.9485 



0.948 



0.9475 



0.9495 



0.949 



0.9485 



0.948 



0.9475 




(21 



X3) 



(a) 



0.10 0.01 3 
0.50 ~ ~ 



0.90 



6 Mq, Z=0.02 
pulse 13 



(b) 



20 40 60 

t (y) to = 74335421.52 y 



80 



Fig. 5. Same as Fig. |, but for pulse 13 of the M6Z02 star 

terpulse phase in order to bring some ^°^Pb to the surface 
of the M3Z02 star. Whether or not such a deep 3DUP can 
operate is still unknown. 

5.2. The M6Z02 case 

In contrast to the M3Z02 case, the M6Z02 star at its 13th 
pulse is the site of the IS+CPS. Neutrons are not produced 
anymore at times t > tin, so that assumption (v) is fully 
valid. 

As in the M3Z02 case, almost no ^"^Pb can be de- 
stroyed during the pulse decay. This results from the fact 
that, at the time of maximum extension of pulse 13, ^°^Pb 
has an effective (mass-averaged) lifetime of the order of 
40 y, while the corresponding ^os^pj lifetime is of the or- 
der of 3 y. These lifetimes have to be compared with the 
approximate 10 y duration of the pulse decay. 
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Fig. 6. Same as Fig. ^, but at time t = 40 y (in the scale shown 
in the abscissa of Fig. ^) following pulse 13 of the M6Z02 star 

The M6Z02 aftcrpulse ^ospb abundance evolution dif- 
fers in two ways from the M3Z02 one. First, the M6Z02 
core mass is higher than the M3Z02 one. This imphes 
much higher intershell temperatures in the former case, 
as demonstrated by comparing Figs. UfaJ andjlfaj or ^(a) 
and Wa). Because of the high sensitivity of A_ to T 
(Fig. g), this leads to ^os^pj ^-decay rates exceeding the 
^°^Pb e~-capture rates over a large fraction of the in- 
tershell region of the M6Z02 star. As a consequence, the 
205p|^ destruction is prevented over more than the inner 
80% of the intershell region [Figs.|(bJ andli'cj]. Second, 
the response time of the envelope to the thermal relax- 
ations following the decline of the pulse is about ten times 
shorter in the M6Z02 than in the M3Z02 case. This favours 
further the survival of ^"^Pb before the eventual occur- 
rence of the 3DUP. 

In conclusion, large amounts of ^'^^Pb can emerge from 
massive AGB stars for moderately deep 3DUP events. 

5.3. The M3Z001 case 

There are many similarities between the 14th pulse- 
interpulse of this star and the M6Z02 case discussed above. 
In both stars, the s-processing is of the IS+CPS type, and 
no neutrons can be produced after the maximum exten- 




-20 20 40 



t (y) to = 279862021.9 y 
Fig. 7. Same as Fig. |, but for pulse 14 of the M3Z001 star 

sion of the pulse. In addition. Figs, ^('bj and ||fcj reveal 
that much of the initial ^°^Pb survives the afterpulse, and 
can be brought to the surface by an eventual 3DUP. 

5.4. The survival o/^^^Pb in the envelopes 

If indeed some ^"'^Pb can find its way to the stellar surface, 
it remains to be seen if it is able to survive there during the 
entire duration tinter of the interpulse phase subsequent 
to the 3DUP. This is a non trivial question as the 205p^ 
e~-capture and ^os^^ /?-decay rates depend on the detailed 
thermal structure of the envelope. 

Table || indicates that no significant ^"^Pb destruction 
occurs in the envelope of the three considered stars, in 
spite of the fact that the effective ^ospj^ decay rate in the 
envelope, obtained by mass-averaging over all its convec- 
tive layers, is increased over its terrestrial value due the 
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0.858 0.86 

Fig. 8. Same as Fig. ^, but at time t = 60 y (in the absissa of 
Fig. |) following pulse 14 of the M3Z001 star 



high temperatures reached at the bottom of the envelope. 
Even so, Table || shows that the effective ^osp]-, lifctii^^^ 
remains much longer than the interpulse duration, so that 
205p|^ has no time to be destroyed significantly in the en- 
velope between two successive pulses. 

6. Discussion 

The following remarks are worth emphasizing: 
(i) the stellar models considered in this work lead to some 
neutron processing through the CPS or IS-I-CPS. Inter- 
pulses that are hot enough for allowing the development 
of a pure IS might well be obtained at later evolutionary 
phases. In such conditions, one might wonder if, follow- 
ing its IS production, ^°^Pb could not be destroyed be- 
fore the episode of maximum extension of the following 
pulse. This would of course invalidate our assumption of 
the presence of the considered radionuclide at this evo- 
lutionary phase and, as a consequence, make the calcu- 
lations reported in the previous sections meaningless. In 
fact, we have performed some numerical tests which sug- 
gest that the IS-produced 2°^Pb, even if transformed into 
^°^T1 after the IS neutron-production episode and before 
the end of the interpulse phase, could be fully replenished 
at an early stage of the following pulse growth by the fast 



^•^^Tl (/3~ ) ^'^^Pb transformation taking place at the high 
temperatures encountered in the bottom layers of the con- 
vective pulse. In such conditions, the situation would be 
essentially the same as the one considered in this work, 
with ^"^^Pb uniformly distributed in the pulse at its max- 
imum extent. Of course, more definite conclusions have 
to await the construction of detailed models leading to a 
pure IS production of ^°^Pb; 

(ii) in the CPS or IS-I-CPS scenarios encountered in this 
work, our ignorance of the precise mass location of the pri- 
mary ^'^C responsible for the neutron release is not likely 
to influence the fraction of ^°^Pb that can survive the ther- 
mal pulses. This results from the fact that the interpulse 
^^C-loaded layers are eflticiently mixed by the time of max- 
imum extent of the next pulse; 

(iii) in contrast, the fraction of surviving ^°^Pb drastically 
depends on the specific conditions found in the intershell 
region during the afterpulse phase prior to the 3DUP; 

(iv) the level of contamination of the stellar surface with 
205pj^ is probably quite insensitive to the precise time 
of occurrence of the 3DUP. Indeed, Figs, ^b), ^(b) and 
^(b) reveal that, in regions where ^"^Pb survives in large 
amounts, it reaches its equilibrium value on timescales of 
the order of the characteristic intershell thermal readjust- 
ment time, i.e. quite likely well before the occurrence of 
the 3DUP; 

(v) the 205p]2 eventually dredged up to the surface would 
not be severely destroyed during the interpulse (Sect. 5.4). 

Of course, the amount of ^°^Pb finding its way to the 
stellar surface does depend on the depth of the 3DUP. 
Figure 9 illustrates this dependence for the M6Z02 model 
displayed in Fig. |[ It is also instructive to present the 
surviving ^°^Pb at the stellar surface versus the surface 
12(-i^16q abundance ratio. Figure 10 shows this correlation 
after a single 3DUP following pulses 12, 13 and 14 in the 
M3Z02, M6Z02 and M3Z001, respectively. As in Fig. 9, 
the depth of the 3DUP is a free parameter. In addition, it 
is assumed that the stellar surface is free of ^"^Pb before 
the considered 3DUP, while the ^^C and ^^O profiles are 
those provided by the models. 



Table 2. Percentage fn'^surv of ^°®Pb in the envelope which 
survives the interpulse of the given star. The is the approxi- 
mate temperature (in 10® K) at the bottom of the envelopes, 
and tinter the duration (in y) of the interpulse. The last col- 
umn gives an estimate of the mass-averaged ^"^Pb e~-capture 
lifetime (in y) in the envelope. 



case 


interpulse 


The 


tinter 


renv 

JPb.surv 'e 


M3Z02 


12-13 


~ 5 


~ 76000 


100 1.2 lO'' 


M6Z02 


12-13 


~ 80 


~ 3000 


99.4 4.8 10^ 


M3Z001 


13-14 


~ 30 


- 10700 


93.3 1.4 lO'^ 
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Fig. 9. Number density of ^"'^Pb and ^^C mass fraction (to 
be read on the left-hand ordinate) and ^'^C/^^O ratio (to be 
read on the right-hand ordinate) at the surface of the M6Z02 
star 40 y after its pulse 13 (see Fig. ^) as a function of the 
parametrized depth of the 3DUP within the ^^C-rich region 
left over by pulse 13. Note that hot bottom burning at the 
base of the envelope is responsible for the low ^^C abundance 
prior to the 3DUP 



Figure 10 demonstrates that substantial amounts of 
205p|-j (2Qui(j pollute the stellar surface well before the con- 
sidered stars become C-type stars, and that the 205p|^ 
amount increases much more dramatically with the in- 
creased depth of the 3DUP than the corresponding ^^C 
surface enrichment. Of course, it has to be stressed that 
the situation displayed in Fig. 10 is quite artificial, as the 
surface ^^C concentration most likely builds up progres- 
sively through successive, and possibly relatively shallow, 
3DUPs. In M ^ 3 Mq AGB stars, in which only very httle 
205p|^ survives in the outer layers of the intershell region, 
the succession of pulses could thus lead to an increase of 
the surface ratio without any associated ^"^Pb 

enrichment. 

7. Conclusions 

This paper presents the first attempt to evaluate in a 
quantitative way the possibility of contamination with 
205p|-, Q-f surface of thermally pulsing intermediate- 
mass AGB stars. We rely on detailed models of three dif- 
ferent stars (M = 3M0,Z = 0.02; M = 6Mq,Z = 0.02; 
M — 3Mq,Z = 0.001), and on some simplifying assump- 
tions concerning in particular the 3DUP characteristics, 
some of which (like its depth) being considered as free 
parameters. From the examination in each of the con- 
sidered stars of a single pulse-interpulse cycle following 
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Fig. 10. Number density of ^'''^Pb as a function of the ^^C/^'^O 
ratio at the surface of the three stars analysed in the text 
following the pulse indicated in Table 1. Each location on the 
curves corresponds to a given depth of the 3DUP, which is 
taken as a free parameter 



a sequence of a dozen of pulses, we conclude that, during 
the afterpulse phase, the ^°^Pb destruction is prevented in 
the hottest (T ^ lO^K) parts of the intershell region. As a 
direct consequence, the chances for a significant ^°^Pb sur- 
face enrichment are likely to increase with M for a given 
Z, and to increase with decreasing Z for a given M. More 
specifically, following the considered pulses at least, the 
enrichment appears to be rather unlikely in the 3 Mq star 
with Z — 0.02, while it seems to be much more proba- 
ble in the other two considered stars. Of course, the older 
the AGB star (i.e. the larger the pulse number), the more 
205pj-, jg likely to survive during the afterpulse, as the tem- 
peratures in the intershell regions increase with the pulse 
number. It may thus be that the surface of a 3 Mq star 
of the M3Z02 type could be polluted with a significant 
^"^Pb amount at a more advanced stage of its AGB evo- 
lution, provided of course that the pulse sequence is not 
too drastically shortened by severe mass losses. 

Quite clearly, the problem of the ^ospb survival in AGB 
stars will have to be revisited more thoroughly once AGB 
models that predict self-consistently the development of 
the s-process and the occurrence of the third dredge-up 
become available. In such conditions, it will be meaningful 
to analyse in detail the surface ^ospi^ enrichment resulting 
from a series of self-consistently computed pulse-interpulse 
cycles. 
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